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We study the anomalous electromagnetic pion form factor _F 7r o 7 * 7 » with a set of holographic 
models. By comparing with the measured value of the linear slope, some of these models can be 
ruled out. From the remaining models we obtain predictions for the low-energy quadratic slope 
parameters of _F^.o 7 * 7 », currently out of experimental reach but testable in the near future. We 
find it particularly useful to encode this low-energy information in a form factor able to satisfy also 
QCD short-distance constraints. We choose the form factor introduced by D'Ambrosio, Isidori and 
Portoles in kaon decays, which has the right short distance for a particular value of the quadratic 
slope, which is later shown to be compatible with our holographic predictions. We then turn to a 
determination of the (dominant) pion exchange diagram in the hadronic light by light scattering 
contribution to the muon anomalous magnetic moment. We quantify the theoretical uncertainty 
in (g — 2) M coming from the different input we use: QCD short distances, experimental input and 
low-energy holographic predictions. We also test the pion-pole approximation. Our final result 
is a I = 6.54(25) ■ 10~ 10 , where the error is driven by the linear slope of i 7 ^o 7 « 7 *, soon to be 
measured with precision at KLOE-2. Our numerical analysis also indicates that large values of 
the magnetic susceptibility xo are disfavored, therefore pointing at a mild effect from the pion off- 
shellness. However, in the absence of stronger bounds on \o, an additional (10 — 15)% systematic 



uncertainty on the previous value for a M cannot be excluded. 



PACS numbers: 14.60.Ef, 13. 40. Em, 11.25.Tq 
I. INTRODUCTION 

The anomalous magnetic moment of the muon 
is presently one of the most stringent tests of the 
standard model and, due to the current experimen- 
tal precision, a very useful tool in the search for 
new physics. As of today, a discrepancy of roughly 
three standard deviations persists between the stan- 
dard model estimate and the experimental value at 
Brookhaven pj. Whether this is to be ascribed to 
evidence for new physics depends crucially on the 
reliability of the standard model computation. 

With the electroweak corrections under good con- 
trol, the main source of uncertainty on the theo- 
retical side comes from the hadronic contributions, 
namely the vacuum and the light-by-light contri- 
butions. While information on the hadronic vac- 
uum polarization can be extracted from experi- 
mental data on the hadronic e + e~ cross section, 
the hadronic light-by-light (HLBL) contribution can 
only be estimated through nonperturbative theoret- 
ical methods, i.e. lattice simulations and hadronic 
models. Due to the complexity of the calculation, a 
lattice result is unfortunately still not available and 
one has to resort to the latter. 

In the past there have been plenty of determina- 
tions based on particular models like the ENJL Q, 
HLS Q and more recently on different models in- 



spired by the 1/N C expansion under the assumption 
of lowest meson dominance (LMD) [l-Q . Different 
as they are in nature, all the abovementioned models 
eventually determine their free parameters impos- 
ing constraints, either theoretical or experimental. 
Therefore, the source of uncertainty generally de- 
pends on: (i) the experimental input used; (ii) the 
different theoretical constraints imposed; and (iii) 
an intrinsic model dependence, which is very dif- 
ficult to estimate. While it is certainly reassuring 
that, ever since the seminal paper of Ref. Q settled 
an important sign issue, the results for the different 
models lie in the same ballpark and no big disagree- 
ment is present, there are certain unresolved issues 
that might affect the theoretical error estimate (see 
for instance the recent discussion in J8j). With the 
proposals at Fermilab Q, J- PARC and poten- 
tially at Frascati (Tlj . the increase in experimental 
precision will aim at reaching 1.5 ■ 10~ 10 accuracy. 
Therefore, a reliable estimate of the theoretical er- 
ror in the hadronic light-by-light contribution be- 
comes of the utmost importance (see, for instance, 
the recent assessments of the error budget given in 
Refs. [Hill). 

In this paper we propose to use the holographic 
principle [14J to address such issues. Holographic 
models of QCD, while also based on the 1/N C expan- 
sion, offer several advantages over four-dimensional 
LMD models. First and foremost, AdS/QCD mod- 
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els are implemented at the Lagrangian level, and 
therefore Ward identities between correlators are au- 
tomatically fulfilled. Second, due to the AdS met- 
ric, short-distance matching to asymptotically free 
QCD is easily achieved. Third, hadronic resonances 
in AdS/QCD models arise as the Kaluza-Klein states 
in the process of compactification from five to four 
dimensions. Therefore, as opposed to LMD mod- 
els, the full tower of states is automatically imple- 
mented and its separate contributions can be ana- 
lyzed. Finally, even though an infinite number of 
resonances is present, the number of free parame- 
ters is very small. Thus, unlike phenomenological 
hadronic models, very little input is needed in order 
to be predictive. Holography, therefore, has the po- 
tential to become a consistent hadronic model for all 
the hadronic light-by-light contributions to [g — 2)^, 
at least at leading order in the 1/N C expansion, i.e., 
in the limit of single-resonance exchange. 

In this paper we will concentrate on the neutral 
pion exchange contribution. Traditionally, this con- 
tribution has been extracted from the knowledge on 
the electromagnetic pion form factor F 7r o 7 . 7 *. The 
questions we want to address are the following: 

(i) Which are the parameters entering i 7 ' 7r o 7 * 7 * 
that mostly affect the uncertainty on (g — 2)^? 

(ii) It has been recently pointed out [l5| that the 
quadratic slope $, defined by 
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might be of importance to reduce the uncer- 
tainty in the (g — 2) M . However, an experi- 
mental determination thereof is still lacking. 
Can holographic models of QCD be used to 
estimate the linear and quadratic slopes of 
i 7 ' 7I .o 7 » 7 *? Do these models satisfy the required 
short-distance constraints? How big is the im- 
pact of the holographic slopes on (g — 2) M ? 

(iii) Is vector-meson dominance a justified approxi- 
mation for i 7 ! lr o 7 » 7 «? How does this approxima- 
tion carry over to the evaluation of the (g— 2) M ? 

(iv) What is the impact of not restricting the neu- 
tral pion to be on-shcll in the HLBL contribu- 
tion? In other words, what is the accuracy of 
the pion-pole approximation? 

In order to address these questions we will study 
_F! n .o 7 . 7 « in a set of holographic models, which mostly 
differ in the way they implement chiral symmetry 
breaking and the spacing in the hadronic vector 



spectrum (Regge-like or not). We will test the mod- 
els against existing experimental data on F 7r a 1 * 1 *- 
In particular, we will compare the different predic- 
tions for the linear slope a with the results from 
CELLO [l6|]. Agreement with experiment will be 
used as a filtering criteria for the different models. 
Then, from the accepted models, we will extract a 
prediction for the quadratic slopes (3 and 7. 

For the analysis of the (g — 2)^ we will adopt a 
strategy similar in spirit to Ref. [l7j], where a set 
of simple interpolators were tested against experi- 
ment and then used to estimate the contribution to 
the HLBL. A particularly useful form factor, which 
displays immediately the required short and long- 
distance properties is found to be the one used by 
DAmbrosio, Isidori and Portoles (DIP) in kaon de- 
cays [lj]: 
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In particular, to study the pion off-shellness we will 
promote the DIP form factor to an interpolator valid 
for arbitrary pion momentum. This leads to a new 
short-distance constraint which is naturally im- 
plemented in our DIP interpolator. 

This paper is organized as follows: in Section [TT] 
we will briefly review the holographic principle and 
its most common realizations to study QCD. These 
models will then be used in Section Hill to study the 
7r°7*7* form factor, with special emphasis on its low- 
energy predictions. In Section llVI we will introduce 
the DIP interpolator and i) give a first estimate of 
the pion exchange contribution to the HLBL piece of 
(g — 2)n and then ii) present an extension of the DIP 
interpolator with an extra pole, which will allow us 
to play with the whole set of long and short-distance 
constraints and test the stability of our results. Fi- 
nally, conclusions will be given in Section [V] Tech- 
nical details are provided in two Appendices. The 
paper is organized such that readers interested in 
phenomenological applications can skip Section [TT] 
without loss of continuity. 



II. HOLOGRAPHIC MODELS OF QCD 

The AdS/CFT conjecture Q offers one of the 
most promising ways to study gauge theories in 
strongly coupled regimes through their weakly cou- 
pled supcrgravity duals compactificd in AdS.g. In 
the recent years attempts have been directed to- 
wards using the gauge/string duality to QCD, i.e., 
a holographic equivalence is conjectured between 
four-dimensional strongly coupled QCD at large N c 
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and a five-dimensional weakly interacting gauge the- 
ory coupled to gravity on a five-dimensional space 
not necessarily (but asymptotically) AdS5. The 
fact that QCD is conformally invariant in the deep 
Euclidean makes the gauge/string duality a good 
starting point towards a theory of hadrons. How- 
ever, QCD is not conformal in the strongly cou- 
pled regime and crucial ingredients like confinement 
or chiral symmetry breaking have to be incorpo- 
rated. Confinement can be easily modeled by mak- 
ing the bulk space compact, for instance by placing 
an infrared brane some distance zq away from the 
ultraviolet brane. The QCD resonances are then 
the Kaluza-Klein modes arising from the compact- 
ification. Chiral symmetry breaking is more in- 
volved but AdS/CFT seems to have the potential 
to describe both explicit and spontaneous symme- 
try breaking The pion field can also be in- 
corporated, and the agreement of the whole picture 
with QCD, especially with vector mesons, is quite 
remarkable |20l - l23l ]. 

In order to bridge the gap between the original 
AdS/CFT conjecture and AdS/QCD two main ap- 
proaches have been followed, the so-called top-down 
and bottom-up. In the top-down approach, one 
looks for a suitable setting of D-branes in string 
theory, with gauge theory on their world volume, 
which at low energy would produce an effective back- 
ground geometry for the dual five-dimensional grav- 
itational theory. An example is the Sakai-Sugimoto 
model [22| to be considered later. 

The bottom-up approach is more phenomenolog- 
ically oriented and one starts directly from warped 
five-dimensional models, with an AdSs metric in the 
ultraviolet regime but with drastic deviations in the 
infrared, where nonperturbative effects of QCD force 
a description in terms of new low-energy degrees of 
freedom. In this paper we will consider Hard- Wall 
(HW) and Soft-WalHSW) models. Both in the HW 
models of [2(| and 24 1, the AdSs space is cut off at a 
finite size, producing an infinite number of Kaluza- 
Klein resonances, to be identified with the hadronic 



spectrum. The two models differ in the implementa- 
tion of xSB, though. In (20|, xSB is induced by a 5D 
scalar field, holographically dual of the qq operator 
of QCD, whose nonvanishing vacuum ex pec tation 
value is responsible for xSB. In contrast, in [2J], xSB 
is achieved by imposing appropriate infrared bound- 
ary conditions. The SW model proposed in [25[ is 
a five-dimensional holographic model in which the 
AdSs space is noncompact. Confinement and an in- 
finite number of bound states follow from the pres- 
ence of a nontrivial dilaton background. The main 
feature of the SW model is its ability to produce an 
hadronic spectrum with Rcggc behavior, which leads 
to better agreement with the resonance spectrum in 
QCD. 

In any of the models above, one boldly conjectures 
the applicability of the holographic recipe to com- 
pute correlation functions of the dual 4-dimensional 
theory pril[27j. For every quantum operator O(x) in 
QCD, there exists a corresponding bulk field <p(x, z), 
whose value on the ultraviolet brane, <j)(x, 0) = 
4>o{x), is identified with the four-dimensional source 
of 0{x). Hence, the generating functional of the 
four-dimensional theory can be computed from the 
five-dimensional action evaluated on-shell (neglect- 
ing stringy corrections), i.e., 



exp (iS 5 [c))o(x)}) 



(cxp 



i / d 4 xcj) (x)O(x) 
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(1) 

Integrating by parts the quadratic part of the five- 
dimensional action on-shell effectively reduces to a 
boundary four-dimensional term quadratic in 4>o(%)- 
By varying the action with respect to 4>o(x) one 
can generate the different connected n-point Green's 
functions of QCD. 

The simplest five-dimensional action can be gener- 
ically written in the form 



5s — 5ym + Sx + Scs , 



(2) 



where 



Sym 



tr / d 4 x I dz e -*«^-3u;(jz) 



^(L) J~{L)MN + J~{R) -r(R)MN 



-pMN 77 



Sx = tr J d A x ° dz e-* {z) w(zf[D M XD M X^ + V(X^X)] 



(3) 



with Tmn = 9mBn— 9nBm— i[&M, Bn] and£>L,# = vector) five-dimensional fields. If these are coupled 
V T A, where V(A) <S U{2)v(A) are vector (axial- to Dirac currents, then the holographic prescription 
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predicts that they are massless and hence gauge- 
invariant. It is common to work in the axial gauge 
V5 = A$ = 0. The surviving five-dimensional gauge 
fields Vfj,(x,z) and A fl (x,z) holographically corre- 
spond to four-dimensional vector and axial-vector 
QCD currents, qj^q and 9^7^75(7^; respectively. 

$(z) is a dilaton field and X(x, z) a scalar field 
transforming under the chiral group as gi,Xg R . 
Then, accordingly, DmX = 8mX — iLmX + iXRm- 
As we shall see below, the presence of the scalar field 
X and the form of its potential depend on the model 
considered. For instance, in models where x^B is in- 
duced by chirally-asymmetric boundary conditions 
on the gauge fields, X(x, z) is not essential. 

The extra-dimension is taken to extend over the 
interval (0, Zo), where the upper limit can be infinite 
for some models. The metric of the five-dimensional 
space can be written generically in terms of a warp 
factor w{z) as 

g M Ndx M dx N = w(z) 2 (r^dx^dx 1 ' - dz 2 ) , (4) 

where X]^ = Diag (1, -1, -1, -1), fj,,u = (0,1,2,3) 
and M 7 N = (0, 1, 2, 3, z). In AdSs space the warp- 
ing factor takes the form w(z) = 1/z. 

As pointed out in (28[, anomalous processes in 
four dimensions can be reproduced from the five- 
dimensional Chern-Simons term 

(5) 

In order to account for chirality, one should work 
with 

S£I S [B L , B R ] = S cs [B L ] - S CS [Br] . (6) 

In the following we will briefly review the general 
features of the holographic models we will consider 
for our analysis of the electromagnetic pion form fac- 
tor, mainly to fix our notation. Further details can 
be found in the original literature. 



A. Hard- Wall models 

The distinguishing features of HW models are a 
compact fifth dimension, < z < zq and a constant 
dilaton field $(z). We will concentrate on two such 
models, which we name HW1 [H[H| and HW2 [H|. 
The presence of the upper bound z = zq defines an 
infrared brane, producing an explicit breaking of the 
scale invariancc of the dual four-dimensional theories 
at energies « 1/zo w 1 GcV. Moreover, in the pres- 
ence of this cutoff, Wilson loops follow an area-law 
behavior in the infrared regime, thereby simulating 
the onset of a confining phase. 



In compliance with x>3B in the four-dimensional 
theory, axial gauge invariance is broken in the five- 
dimensional background, the longitudinal compo- 
nent of the axial- vector field becoming physical and 
related to the pion field. 

HW1 and HW2 differ in the mechanism leading 
to the xSB. In the HW1 model, the breaking is due 
to the scalar field X(x,z), whose coupling to the 
axial-vector gauge field produces an effective mass 
term, breaking the gauge invariance in the axial sec- 
tor. In the HW2 model, instead, no such field is 
present in the five-dimensional Lagrangian @ and 
%SB is achieved through different infrared boundary 
conditions for vector and axial- vector fields. 



1. x^B from a scalar bulk field: HW1 

One considers the action Q without dilaton 
field and with the complex scalar field X(x, z) = 
v(z)U(x, z)/2, where U contains the pion field and 
v(z) is the scalar component that breaks chiral sym- 
metry in the bulk. The scalar potential is reduced 
to the mass term, given by: 

V(X^X) = ^X*X . (7) 

The previous expression guarantees that X can be 
coupled to qq. Solving the five-dimensional equa- 
tions of motion at zero four-dimensional momentum 
one finds v(z) = (m q z + crz 3 ), where the parame- 
ters m q and a are holographically identified as the 
quark mass and the qq condensate, i.e. the sources 
of explicit and spontaneous x^B, respectively. 

A nonvanishing v(z) induces a z-dependent mass 
term for the axial- vector field in Eq. (j3J), thereby 
breaking the degeneracy between vector and axial- 
vector resonances. The scalar field X also pro- 
duces a nontrivial coupling between the longitu- 
dinal component of the axial-vector field and the 
pion field. Defining A^Jq, z) = —iq ll Lp{q,z) and 
U(x,z) = exp[2it a ir a (x,z)} (for SU(2), t a = a a /2, 
with a a being the Pauli matrices) one gets the sys- 
tem of coupled equations [2(| : 

dJ- z d z ^+^^(^-^) =0, (8) 

^ dzlp a + 9_t^l dz7r a = (9) 
Z Z 

As a result, the HW1 model satisfies the Gcll- 
Mann-Oakes-Renner (GMOR) relation [29j : 

mlf 2 = (m u + m d ) < qq >= 2m q a . (10) 

The pion wave function ip(z) = 1 — ^(z) can be ob- 
tained from the coupled equations ([5]) and © which 
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describe the dynamics in the axial sector. In the 



chiral limit, 



= 0, one obtains 



(ii) 



where £ = a being the quark condensate. 

The pion wave function is different from zero on the 
infrared branc, 



*(zo) 



1 



V3r[2/3] 
7rI- 2 /3(^o 3 ) V2(^o 3 ) 2 



1/3 



(12) 



a fact that will be important in the evaluation of the 
anomalous amplitude. 

By solving the equation of motion for the vector 
field (Q 2 = —p 2 being the 4-dimensional Euclidean 
momentum): 



z I z 







(13) 



subject to the (Dirichlct) ultraviolet boundary con- 
dition S(Q, 0) = 1 and the (Neumann) infrared one, 
d y J(Q, zq) — 0, one finds the so-called vector bulk- 
to-boundary propagator S(Q, z), which can be writ- 
ten in terms of Bessel functions: 



J(Q,z) = Qz 



K 1 {Qz)+h{Qz) 



Kq(Qz ) 
Iq(Qzo) 



(14) 



Vector resonances are associated with solutions 
tpn(z) of the equation of motion for the vector field 
at discrete p 2 - 
ditions: tp n (0) 



m„, with vanishing boundary con- 
0, and d z tp n (z ) = 0. 



2. X^B through boundary conditions: HW2 

The action of the HW2 model [24| is entirely given 
by the gauge field part of eq.©. The role of the 
scalar field X as the source of chiral symmetry break- 
ing is played by asymmetric boundary conditions be- 
tween vector and axial fields. In this model, vectors 
are required to obey infrared Neumann boundary 
conditions while axial fields satisfy Dirichlct bound- 
ary conditions. This leads to a splitting in the mass 
spectra that qualitatively reproduces what is ob- 
served in nature. 

The pion field is built from Wilson lines extending 
between the 5D boundaries: 



U(x)=U(x)^ 1 (x) , (15) 
Zl,r( x ) = Pexp|-z^ ° dz'Bi' R (x,z')\ . (16) 



where 



The fact that the pion is a nonlocal object leads to 
difficulties to implement the GMOR relation. For 
our purposes, this will not be of relevance, since we 
shall consider the HW2 model only in the chiral limit 
of massless pions. 

Chiral symmetry breaking is implemented by 
splitting the fields like 

% (x, z) = Vfj, {x, z) + % (x, 0) , (17) 
(aj, z) = (x, z) +a (z) (x, 0) , 

where the last terms are the sources and cx(z), which 
plays the role of the pion wave function, is deter- 
mined by demanding no mixing between the pion 
and the axial resonances: 



a(z) = I--, 



(18) 



For comparison, we will also consider the flat case, 
i.e. w(z) = 1. In that case the bulk-to-boundary 
propagator takes the form 

J(Q, z) = cosh(Qz) + tanh(Qzo) sinh(Qz) , (19) 

and the pion wave function is given by 



a{z) = 1 



ZQ 



(20) 



3. The Sakai-Sugimoto as an HW2 model 

The original action of |22j is 

S = Sym + Scs j (21) 



with 



dz tr 



(22) 

where the functions h{z) and k(z) are given by 

h(z) = (1 + z 2 )- 1 / 3 ; k(z) = l + z 2 , (23) 

and 5cs is given in Eq. ([6]). 

The constant k is related to the 't Hooft coupling 
A and the number of colors N r as 



(24) 
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The model also has a mass scale Mkk which in (|22|) 
was absorbed in the dimensionless parameter z. For 
our numerical analysis, the two parameters will be 
chosen as 122l 



M K k = 949 MeV , 



0.00745 



(25) 



G 



in order to fit the experimental values of the p meson 
mass, m p ~ 776 MeV, and the pion decay constant, 
~ 92.4 MeV. 

The action (j2"Tj) was obtained in [22[ as the ef- 
fective action of Nf probe D8-branes placed in 
the background of N c D4-branes, studied in (30j . 
and stands for an effective theory of mesons in 
four-dimensional (large N c ) QCD with Nf mass- 
less quarks. In the following we will reformulate the 
Sakai-Sugimoto model and show that it can be cast 
in the form of a HW2 model. This reformulation 
will prove useful for computational purposes in the 
following Section. 

As a first step, let us define y = tan -1 z, with 
values in the finite interval — tt/2 < y < ir/2, whose 
end points y = ±tt/2 correspond to the ultraviolet 
branes while the point y = corresponds to the 
infrared brane. In the A z = gauge, (|2"2"j) becomes 



'YM 



d A x 



r/2 



dy tr 



-tt/2 



h(y)Tl-Tl 



(26) 

where h(y) is an even function of y. This allows us 
to decompose the gauge field in parity-even and 
parity-odd parts as Ap = V^+A^, with V^(x 7 —y) = 
Vfj,(x,y) and A M (x, —y) = -A^(x,y), and restrict 
the theory to half the interval, i.e. < y < tt/2. 
Parity transformations in the fifth dimension corre- 
spond to the exchange of left and right-handed chi- 
ralities and can be used to distinguish between vec- 
tors and axial- vector fields (vectors and axial- vectors 
having even and odd (/-profiles, respectively). As a 
consequence, they satisfy different boundary condi- 
tions on the infrared brane (at y = 0), i.e. Neumann 
for vectors, dyV^y^ = 0, and Dirichlet A„\, l= n = 



n\y=o 



for axial-vectors. 

Next, let us define the dimensionless variable z = 
tt/2 — y, with < z < tt/2 (not to be confused with 
the original variable in (|2"21 ). z = corresponds to 
the ultraviolet brane and z = ±tt/2 to the infrared 
brane, and the action takes the following form: 

tt/2 

dztr|(sinz)- 4/3 

x [(V/ M „ - i[Af+, A„]) 2 + {DpAv - DvArf] 
-2(d z V„) 2 -2(d z A^} , (27) 

where V^ v = d^V v - d v V^ - V v ] and D^A V = 
duAv-ilV^Av]. 

In the new variable z <G [0, tt/2] the SS model takes 
the form of the HW2 model with different boundary 
conditions for vector and axial-vector fields on the 
infrared boundary z = tt/2, and an effective metric 
which is not AdSs: 

gMN<ix M dx N = (smz)~~5r] tl vdx^dx L ' — (sin z)~ ^dz 2 . 

(28) 



The chiral field containing the pion is also built from 
Wilson lines extending between boundaries. One 
can easily show that in the SS model the function 
a(z) discussed in the previous subsection takes the 
simple form 



a(z) = 1 



2z 

IT 



(29) 



Vector (and axial-vector) fields can be obtained by 
solving the bulk-to-boundary differential equation 



Q' 2 



dz 2 ' (sinz) 4 / 3 



J = , 



(30) 



subject to the boundary conditions Jv,a(QiO) = 1, 
and d z J v (Q,TT/2) = 0, J A (Q,n/2) =0, for which 
there is no analytic expression. In Fig. 1 we show 
the first three normalized eigenfunctions for vector 
and axial- vector fields in the new extra-dimensional 
variable z. As usual, the Fourier transform of 



~ u t 

N , V 

I I 
\ \ 

-0.2 * * / 
\ \ / 
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FIG. 1: The first three vector and axial-vector meson 
resonances of the SS model, using the new variable z £ 
[0,7r/2]. 

the five-dimensional gauge field has been written 
as A^(q,z) = A^,(q)ip(z), with ip(z) satisfying the 
equation 



T 



(sinz) 4 / : 



:1)>{Z) = . 



(31) 



Dynamical resonances ip n (z) then correspond to nor- 
malized solutions of ([3"1"T) with q 2 = m^, where the 
normalization condition is given by 

,7r/2 

4 K / (sinz)- 4 / 3 V> m (z)V>„(z) = S mn , (32) 
Jo 

to ensure the canonical four-dimensional kinetic 
term for A^x); the discrete spectrum of eigenval- 
ues q n = m n , with the corresponding eigenfunctions, 
can then be obtained numerically. 
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B. Soft- Wall models 

SW models were originally motivated as holo- 
graphic models with vector resonances displaying 
Regge trajectories. The action of SW models is given 
by where the metric is AdSs, and e - *' 2 * 1 rep- 
resents a nontrivial background dilaton field given 
by [H 



2 2 
K Z 



leading to a spectrum given by 

m 2 = 4K 2 (n + 1) . 



(33) 



(34) 



The dimensionful constant k can be fixed by fitting 
the mass of the first vector resonance to that of the 
p meson, i.e. k = m p /2. Therefore, the equations 
of motion for the different fields can be obtained 
from those of the HW1 model by replacing the AdSs 
warp factor with e~^^/z. Contrary to the HW 
models, in the SW model the extra dimension is no 
longer restricted to a finite interval, i.e. zq = oo 
and the infrared boundary conditions are replaced 
by normalization conditions on wave functions and 
bulk-to-boundary propagators. 

For the vector fields one has to solve the following 
equation: 



d z j) -Q 2 



-J = 



(35) 



whose solution can be cast in the integral represen- 
tation 13111 



J(Q,z) = k 2 z< 



(l-xf 



exp 



1-x 



2 2 
K Z 



dx , 
(36) 



where a = Q 2 /4k 2 . 

One of the main drawbacks of the orginal SW is 
that chiral symmetry breaking is not implemented in 
a satisfactory way. Due to the absence of an infrared 
brane, the parameters for explicit and spontaneous 
symmetry breaking are not independent. This not 
only invalidates general relations like the GMOR, 
but makes the whole pion dynamics unclear. Recent 
proposals have tried to circumvent this problem, so 
far only at an heuristic level. In this paper we will 
adopt the prescription given in [321 ] . where the pion 
wave function is assumed to be Gaussian, 



a(z) 



(37) 



A comparison of the pion wave function profiles for 
the different models is shown in Fig. 2. 
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FIG. 2: Comparison of the pion wave functions of the 
HW1 model (red solid line), the HW2 model (black dot- 
dashed line) and the Gaussian ansatz of Eq. (|37[) (blue 
dashed line). The vertical line corresponds to the in- 
frared brane at z — zo- 



III. 



7T°7*7* 



FORM FACTOR 



In this Section we will apply the methods de- 
scribed previously for the anomalous pion form fac- 



tor. We will define the 



7r°7""7' 



form factor as 



d'x e^*(0\T{Jg M (x) J£ M (0)} \*°(p)) (38) 



where p = q\ + (72 is the pion momentum, q\ , 52 are 
the momenta of photons, and q\ 2 = — Q\ 2 . Notice 
that our conventions differ by a sign from the ones 
used in [alEl. 

Even though the form factor cannot be computed 
from first principles, there are certain kinematical 
limits where theoretical or experimental information 
is available. For instance, when both photons are 
on-shell, the form factor is determined (in the chiral 
limit) solely by the WZW anomaly term 



7*7* 



,0(0,0) = - 



N c 



12tt 2 U 

and it is therefore convenient to define 

N c 



(39) 



F 7 » 7 * 7r o (Q 2 , Q 



12n 2 f n 



K(Q 2 ,Q 2 ), (40) 



such that K(0, 0) = 1. 

On the experimental side, studies have focused on 
the low-energy behavior, when one of the photons 
is exactly on-shell and the other slightly off-shell. 
In this kinematic regime it is common to define the 
slope of the anomalous form factor a n as 



K(0,Q 2 



a7r -Q 2 
mi 



7T 



. (41) 
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i.e., 



dA(Q 2 ,0) 



dQ 2 



(42) 



Q 2 =0 



The world average is presently = 0.032(4), mainly 
driven by the results of the CELLO Collabora- 
tion [ll| on 7T° — > e + e~7. Given the importance 
of this kinematical regime for the evaluation of the 
(g — 2)^, we also define the curvature as 



dK(Q 2 ,0) 



dQ 4 



(43) 



2 =o 



Later on we will use the different holographic models 
to obtain a determination for both a T and b n . 

Another source of information comes when one of 
the photons is on-shell and the other far off-shell. 
Then the expected behavior of the pion form factor 
is dictated by the Brodsky-Lepage quark-counting 
rules 13a [3611 . 



1 



lim K(0,Q')~ — 



(44) 



Data on this kinematical regime is available from 
e + e~ — > e + e~7r°. The behavior of the previous 
equation is reasonably reproduced by CLEO 13711 
but incompatible with recent data by BABAR |38l |. 
Hopefully upcoming BABAR data with better 
statistics will help clarify the situation. While the 
discrepancy is surprising, we want to emphasize 
that the kinematical regime we are considering does 
not accept an OPE expansion and, therefore, that 
Eq. (031 is not on the same footing as a short- 
distance constraint. 

Short-distance constraints can be obtained when 
both photons have large and equal virtualities. In 
this case, 



lim K(Q 2 ,Q 2 ) 



N c Q 2 



(45) 



It is worth stressing that all the constraints and 
experimental information presented above refer to 
the pion being strictly on-shell. In our case, since 
we are working in the chiral limit, to a very good 
approximation p 2 = m 2 ~ 0. We will come back 
to the on-shellness of the pion when we discuss the 
pion exchange diagram in (g — 2)^. For the time 
being, we will concentrate on the pion form factor 
as predicted from the different holographic models. 



A. Holographic predictions 



In terms of holographic QCD, K(Q\,Q\) can be 
obtained from the VVA terms of the Chern-Simons 



action. Working in the axial gauge, B z = 0, and 
gathering the relevant pieces trilinear in the fields 
one finds: 



N f 
—Eje^tr / d A xdz{d z B^) 

Z47T J 



-<AdS 

'cs 

(46) 

Replacing A p = 9 M 7r and taking the pion to be on- 
shell one ends up with 



S AdS = J*!L. e l»>P° f d 4 x f° dz 



12tH 



(47) 



d z p d p v; d a % - f3d z (d p v« d a v v 



where /3(z) stands for the pion wave function, which 
we will denote as ^f(z) or a(z), depending on the 
model. 

Integrating by parts over z in the second term 
above and dismissing a boundary term (cf. [13] ) in 
order to reproduce the standard four-dimensional 
WZW action one gets 



Nr. 



S£f = 4^ P 7 n dz (d z (3) I d*x n a (d p V«) [d a V v 

(48) 

Variation of Sq^ s above gives the three-point 
function: 



Nc Pa p a tz I rvi r\1\ 



12tt 2 p 2 



(49) 

In compliance with the large-A'c limit, in extra- 
dimensional models the form of the interaction ver- 
tex between the pion and two external electro- 
magnetic currents is mediated by the full tower of 
Kaluza-Klein vector mesons. From the variation of 
the Chern-Simons action with respect to the exter- 
nal sources one obtains 



K(Qi,Q z 2 ) = J(Q 1 ,z)J(Q 2 ,z)d z (3(z)dz 
Jo 

+ (possible boundary terms) , (50) 

where ^7(<5, z) is the vector bulk-to-boundary prop- 
agator at Euclidean momentum Q 2 = —q 2 , defined 
in the previous Section. 

The significance of boundary terms in the Chern- 
Simons action was studied in detail in (28|. In the 
context of the HW1 model, boundary terms in the 
expression of K{Q\,Q\) were shown to be needed 
in order to have the right normalization required by 
the QCD axial anomaly, namely A"(0,0) = 1. The 
presence of these boundary terms is closely related 
to the infrared behavior of the pion wave function. 
As pointed out in Eq. (TT2|) . in the HW1 model the 
value of the pion wave function at zq does not cancel 
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and requires the addition of a boundary term [3J: 



J(Qi,z)J(Q 2 ,z)d z ^(z)dz 



+J(Qx,z )J(Q 2 ,zo)*(zo) ■ (51) 



Indeed, as pointed out in [34| . a(z) and the pion 
wave function (jlip coincide in the deep ultraviolet, 
since at small z 



d z ^(z) ~ -f 2 g 2 z 



d z a(z) 



(54) 



No boundary term is needed both in the case of 
the HW2 and the SS models. The correct normal- 
ization is obtained in both cases (using J7(0, z) = 1) 
due to the boundary conditions satisfied by a(z): 



#(0,0) 



d z a(z)dz = a(0) = 1 



(52) 



For the SW models, since the pion wave func- 
tion is introduced by hand, it depends on the ansatz 
made. In general, pion wave function that do not 
cancel at the infrared brane will require boundary 
terms to correctly implement the axial anomaly in 
K(Q X ,Q 2 ). 

Equipped with Eqs. ([50|) and ([ST]) , together with 
the different expressions for the pion wave function 
and vector bulk-to-boundary propagator, one can 
easily find the expressions for the pion form fac- 
tor. We are mostly interested in checking whether 
short-distance constraints are satisfied and to what 
extend the low-energy information complies with ex- 
perimental data. 



B. Large-Q 2 behavior 



Notice that differences between HW1 and HW2 only 
affect the infrared of the theory, and therefore are 
of no relevance for the large-Q 2 behavior (at most 
exponentially suppressed like 0(e~Q z °)). 

The form of the coefficient in front of Eq. (|53j) 
will obviously change depending on the details of the 
model. For instance, the SW model with the ad hoc 
Gaussian pion wave function is recovered by replac- 
ing Zq — > 1/k. Since zq is related to the pion decay 
constant, while k is matched to the p(770) mass, if 
we impose numerical agreement between both mod- 
els we find the relation 



K = 8tt 2 / 2 



(55) 



which can be compared for instance with the predic- 
tion m 



16^6 2 2 



(56) 



coming from a large- N c sum rule analysis of Tlvv 
and 11^. It is reassuring that both predictions are 
in excellent agreement. 

When Qi = Q 2 = Q, Eq. ^ simplifies to 



In asymptotically AdS holographic models, it is 
easy to show that the function K(Qi, Q 2 ) automat- 
ically satisfies the high-energy constraints discussed 
in the previous Section. For instance, for arbitrarily 
large Qi and Q 2 , and working in the HW2 model, 
it is not difficult to show that the general expression 
reads 



K(w,Q 2 ) 



z 2 Q 2 



x / dtid (VT+w£) k x (VT=^g) £ 3 

Jo 

2 

- ,„3 r 2 n2 [^-(l-w 2 )tanh- w] , 



(53) 



where w = 



2Q 



& 2Q 2 = Qi 



Q\ and £ = 



It is easy to show that the momentum dependence 
of the previous result also holds for the HW1 and 
SW models: AdS dictates the large-Q 2 behavior of 
the bulk-to-boundary propagator, and this is com- 
mon to all of them. Moreover, the shape of the pion 
wave function is very similar for all the models close 
to the ultraviolet boundary, as illustrated in Fig. 2. 



K(Q 2 ,Q 2 



3^2 



(57) 



With the expression for zq in the HW2 model, z$ = 
2/g 2 f n , and g 2 = 12n 2 /N c , it is easy to show that 
it agrees with Eq. (|45|) . Actually, one can go even 
further to show that the leading terms for A — > oo 
are 



K{\ 2 Q 2 ,{\ 2 Q 2 - P 2 )) ^ 

2ff 8 VM 1 , 1 P 



3 Q 2 I A 2 A 3 Q 2 



(58) 



which, up to 0(as) corrections which cannot be cap- 
tured by the HW2 model, is the short-distance be- 
havior found in 0] and (40j . 

One can also explore the regime where one photon 
is on-shcll and the other far off-shell. In that case, 



K(0,Q 2 ) 



•2n2 



zlQ 2 



Sn 2 f 2 
Q 2 ' 
(59) 

which displays the Brodsky-Lcpagc scaling, cf. 
Eq. dll|. 
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Eqs. H53J), ([57]), {55J and ([5ji do not hold, how- 
ever, for the models without asymptotic AdS metric, 
i.e the Sakai-Sugimoto and the flat HW2 model. At 
least for the latter, calculations can be performed 
analytically with the results 



K{w, 



1 



K(Q 2 ,Q 2 
K(0,Q 2 



2wzoQ 

x [vT + w tanh (Qz Vl + w) — {w —> —w}] , 
1 



2Qz a 
1 

Qz 

which fail to reproduce the OPE of QCD. 



(60) 



C. Small-Q behavior and predictions for the 
parameters and b n 



Model 


a (GeV" 2 ) 


(GeV" 4 ) 


7 (GeV" 4 ) 


HW1 


-1.60 


3.01 


2.63 


HW2 (AdS) 


-1.81 


3.65 


3.06 


HW2 (Flat) 


-1.37 


2.25 


2.25 


SS 


-2.04 


4.56 


3.55 


sw 


-1.66 


3.56 


2.76 



TABLE I: Values of a, ft and 7 for the holographic mod- 
els discussed in the main text. 



(hi) SW model: 

z 2 f 1 
g{z) = -— I exp 



,2 f i 



x 



h(z) = 



32k 2 



exp 



1-x 

X 



2 2 
-K Z 



In . 



1-x 



2 2 

K Z 



(1-x) 
In 2 x 



-dx 



■dx 



(65) 



At small virtualities one can expand K(Q 2 , Q 2 ,) in 
the form 

K(Ql Qt)^i + & (Qi + + P Q\Ql 

+7 (Qt + Qt) ■ (61) 

By comparison with the previous Section, we imme- 
diately conclude that = —am 2 and = 7m*. 
The parameter does not contribute to processes 
when one of the photons is real. Even so, it will be 
a useful parameter in the determination of (g — 2) M 
in Section ITVl 

Eq. (|61[) can be reproduced from holographic mod- 
els by working out the small-Q 2 behavior of J(Q, z), 
which we will parametrize as 

Q 2 g(z) + Q A h{z) . 



J(Q,z) 



1 



(62) 



The functions g(z) and h(z) can be easily obtained 
by solving perturbatively in Q 2 the equation of mo- 
tion for the different models. This leads to the fol- 
lowing analytic expressions: 

(i) HW1 and HW2 models: 



h(z) 
(ii) SS model: 
g{z) = M~ 2 K 

h(z) = M- 2 K 



l-21og - 



z 
16 



2(^° 



5 , ( z 
7 + log — 

4 \zq 



(63) 



3 dy (siny)4/3 

, y g(y) 
v (smy)^ 



r/2 



dy 



(siny) 4 / 3 

fV , 2 g(y) 
dy- ' 



(siny) 4 / 3 



(iv) flat case: 

g(z) - 



z 2 


zo 


1" 




z 


~ 2 


z 4 
24 


1 - 


-4^ 

z 



(66) 



Plugging the previous expressions into Eq. (|50[) or 
(|51|) one can obtain the different determinations for 
the low-energy parameters. For the HW1 model one 
finds 



-j-^(z ) + / g(z)d z V(z)dz 
4 Jo 



7 



3z 4 



64 



*(*o) 



h{z)d z ^{z)dz 



16 



<F(z ) - / g{zyd z *(z)dz , (67) 



while for the remaining models: 



P 



7 



g(z)d z a(z)dz , 
h(z)d z a(z)dz , 



g(z) 2 d z a(z)dz 



(64) 



The values for the different models arc collected in 
Table I. Experimentally, only the slope has been 
determined, a = —1.76(22) GeV -2 [4l|, which is 
correctly reproduced by the HW1, HW2 and SW 
models. Notice that, even though the slope is a gen- 
uine low-energy quantity, models without an asymp- 
totic AdS metric fail to reproduce the experimental 
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HW1 1.20 -0.18 -0.04 1.10 -0.06 0.01 1.20 -0.22 0.06 

HW2 (AdS) 1.30 -0.37 0.06 1.10 -0.11 0.01 1.30 -0.37 0.08 

HW2 (flat) 0.99 0.01 0.00 1.00 0.00 0.00 1.00 0.00 0.00 

SS 1.70 -1.10 0.49 1.30 -0.34 0.07 1.60 -1.10 0.54 

SW 0.75 0.14 0.05 0.87 0.09 0.02 0.88 0.09 0.02 

TABLE II: Contribution to a, j3 and 7 due to the first 
three vector-meson radial excitations. For a n , each sub- 
column (n = 1, 2, 3) contains the contribution of k, I < n 
terms. /3„ and 7™ are defined analogously. 




value. The quartic parameters $ and 7, using the 
phenomenologically acceptable HW1, HW2 and SW 
models, are predicted to be 

/3 = 3.33(32) GeV -4 , (69) 
7 = 2.84(21) GeV" 4 . (70) 

Both holographic predictions will be used in the eval- 
uation of the HLBL contribution to the (g — 2)^. 



FIG. 3: Hadronic light by light contribution to (g — 2) M . 

resonances. Finally, let us remark that the repre- 
sentation of Eq. ([72]) is only reliable for moderate 
values of the resonance indices. Most likely Eq. (|7"2"j) 
is, at most, an asymptotic expansion and, therefore, 
beyond a certain threshold it ceases to be meaning- 
ful. 



D. A comment on lowest meson dominance 

Before moving to the evaluation of the (g — 2)^, it 
is important to explore the resonance contributions 
to the results of Table 1. For this, we write the 
vector bulk-to-boundary propagator in its spectral 
decomposition 

^) = Eq4^M, (71) 

ri=l " 

where m n is the mass of the nth vector resonance 
and is related to its decay constant. Plugging 
this expression into Eq. (f5T)f , one obtains the double 
series 

00 

X iix+zm+zis - (72) 

This expression is of the form expected in the large- 
N c limit, with the full tower of vector mesons prop- 
agating between the pion and the two photons. The 
contributions for d, j3 and 7 from the first three 
vector radial excitations are reported in Table II. 
The important point to notice is that lowest vector 
dominance is a feature common to all models (from 
the practically p-dominance of the flat model to the 
more moderate behavior of the SS model). Inciden- 
tally, notice also that, with the exception of the SW 
model, the contribution of the first resonance tends 
to overshoot the total value and has to be compen- 
sated by negative contributions from higher-order 



IV. THE HADRONIC LIGHT BY LIGHT 
CONTRIBUTION TO THE (g - 2) M 

With experimental accuracies at the 10~ 10 level, 
precise determinations of the HLBL contribution 
to the (g — 2) M become of paramount importance. 
The problem is that genuine nonperturbative tech- 
niques are required and having theoretical uncer- 
tainties under control is certainly a challenging task. 
The HLBL contribution to the (g — 2) M is depicted 
in Fig. 3. When internal momenta in the loop 
are high enough one can use perturbation theory, 
but a proper evaluation also requires low energies, 
i.e. Goldstone bosons and hadrons. The neutral 
pion exchange contribution, depicted in Fig 4, turns 

out to be the dominant piece (<4r > ~ 7 ■ 10- 10 ) 

followed by the r\ and rj contributions (a^ l ' r> ' ~ 
3 ■ 1CP 10 ). Quark and Goldstone loops, axial- vectors 
and scalars are also expected to contribute at the 
level of (1— 2)-10 -10 , but cancellations occur (scalars 
and Goldstone loops contribute negatively) and so 
the Goldstone boson exchange ends up collecting the 
bulk of the effect. In particular, this means that the 
naive large— N c counting works and single- resonance 
exchange is the dominant effect. This observation 
suggests that procedures based on the 1/N C frame- 
work are suitable tools to address the problem. 

Different paramctrizations have been used in the 
past to evaluate the HLBL contribution. In par- 
ticular, there is the detailed study of Q with dif- 
ferent paramctrizations based on LMD. Having in 
mind what we did so far, one could be tempted to 
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use the results of the previous Section on the pion 
form factor to evaluate each of the 7r°77 vertices 
in Fig. HI as has been recently done for the HW1 
model [42J. However, in this paper we want to be 
able to test the pion-pole approximation. In order 
to do so, one should trade the ir°jj form factor for 
the more general objects P°jj or A°jj. The issue 
of on-shellness versus off-shellness of thepion has 
attracted some attention in recent years [a 51, 53 . 
Here we will follow the approach taken in [6| and 
characterize the degree of off-shellness entirely by 
the short-distance constraint 



lim F> r7 ,(Q 2 ,Q 2 ,0) 




-Xo 



(73) 



FIG. 4: One of the diagrams contributing to pion ex- 
change in the HLBL contribution to (g — 2) M . 



From the previous expression one readily sees that 
the degree of off-shellness is then regulated by the 
parameter xo- Unfortunately, this parameter is only 
poorly known. We will discuss this issue and its 
impact on the (g — 2) M at the end of this Section. 

Our strategy to test the pion-pole approximation 
will be the following: 

(i) We will use the following ansatz for the form 
factor: 



2 9 



(74) 



first introduced in the context of Kl — > fJ- + fJ-~ 
decays (i~8j . This parametrization was orig- 
inally put forward to study the low-energy 
slope of three-point functions, and has been 
used in experimental studies of Kl 



and Kl — > ir°e + e 



MM 



(ii) We will promote the ansatz in ([71)1 to an an- 
alytic interpolator valid for any value of the 
photon momenta. Then, using Eqs. (HUJ) and 
(|T3")) . we get an expression for xo m terms of 
the linear slope A 



Xo 



4tt 2 /| 



(1 + A) . 



(75) 



Aside from being simple and with a very intuitive 
low-energy behavior, the parametrization of Eq. (|74p 
has additional advantages. One of the most impor- 
tant is that the computation ofjo — 2)^ turns out to 



be greatly simplified. In Ref. [4j it was shown that 
if the interpolator can be cast in the generic form: 



F o . 



Nr. 



12^L l fiql) 

J2^2 — —St(ll) , (76) 



then, in the general two-loop expression for the pion 
contribution to the HLBL scattering, all angular in- 
tegrations can be performed using the hypcrspher- 
ical approach. The diagram of Fig. HJ can then be 
evaluated through the following expression: 



where 



—) {%(!)+ %{2)) . (77) 



V(i) 



and 



V(2) 



dQi / dQ 2 w 1 (Q 1 ,Q 2 )G 1 (Q 1 ,Q 2 ) 
o Jo 



-w 2 (m v , Qi, Q 2 )G 2 {Qi,Q 2 ) 



(78) 



dQi / dQ 2 w 3 (m v ,Qi,Q 2 )G 3 (Qi,Q2) 
o Jo 



+w 3 (m„,Q 1 ,Q 2 )G4,(Q 1 ,Q 2 ) 



(79) 



In the previous expressions, G, are generalized form 
factors and Wi weight factors, whose expressions are 
given in Appendix A. 

In other words, the contribution to the (g — 2) M 
can be reduced to a more tractable double integral. 
It is not difficult to verify that the interpolator of 
Eq. (J71J) is of the form dTSJ) with 



f(l 2 ) 
9v{q 2 ) 



l + A + (A + ?y) — 
q 2 



X + i] 



9 2 

q 2 - m z v 



(80) 
(81) 



Notice that the interpolator of Eq. ([74")) is assumed 
not to depend on the pion momentum. This simpli- 
fying assumption, together with the specific form of 
Eq. (f76|) . is crucial to be able to use Eqs. (|78]l and 

Our parametrization has 3 free parameters, A, i] 
and my which will be determined by combining the 
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following constraints: 



A 



rnr 



2A 
A 



V = 



n 



(-1.76 ±0.22) GeV" 
, 

3m v 



(82) 
(83) 
(84) 



The first one is a low-energy experimental constraint 
on the slope of the pion form factor (the a param- 
eter of the previous Section), while the last two are 
short-distance information, namely the requirement 
that K(Q 2 ,Q 2 ) at high energies does not go like a 
constant (first equation) and that its 1/Q 2 behavior 
has the right coefficient (see Eq. (|4"5j) . with N c = 3). 
Using as input parameters 



m p = 105.658367(4) MeV 

m v = 134.9766(6) MeV , 

U = 92.4 MeV , 

a em = 1/137.03599976 , 



(85) 



the parameters of the model become 

A = -0.73 ±0.05 , 
V = 0.46±S;i° , 



my 
Xo 



(0.64™) GeV, 
(2.42 ±0.17) GeV" 



(86) 



where the uncertainties are due to the experimental 
error on the slope. 

Several comments are relevant at this point: 

(i) The values of $ from the HW1 and HW2 mod- 
els in Table I are in very good agreement with 
the value of r\ (~ $my) in Eqs. ([86|). We have 
already mentioned that also the linear slope, 
A, is very well reproduced by HW1 and HW2 
in Table I. 

(ii) The value we find for the mass scale has to be 
interpreted as an effective measure of the rele- 
vant scale for the problem. Therefore, we con- 
firm that m p is indeed very close to the natural 
scale to estimate the pion exchange contribu- 
tion. 

(iii) The prediction for xo m fl86[) is obtained from 
the requirement that DIP is a good interpo- 
lator also for an off-shell pion and uses only 
the phenomcnological linear slope, A (and not 
the holographic one). However, the result 
would have been very similar by trading the 
phenomcnological linear slope with the holo- 
graphic values predicted by HW1 and HW2. 



(iv) In order to test the stability and improve the 
accuracy of the results on HLBL one can add 
an extra pole to the interpolator. We shall 
follow this line of thought later on. 



With the values given in Eq. (|86p . the prediction 
for the anomalous magnetic moment is 



al =6.7(3) -10" 



(87) 



In Table III the comparison is made with the LMD 
result quoted in Ref. [4j . In LMD models it is com- 
mon to identify my = m p . In contrast, in our case 
the mass scale is determined dynamically by the con- 
straints. The first thing to realize is that the contri- 
butions for each integral (each column in Table III) 
are in good agreement, even though the interpolators 
are very different and my = 0.64 GeV 7^ m p . How- 
ever, if we redo our analysis but now set my = m p 



from the start (and ignore Eq. ([B"2"j)). then there is 
full agreement with LMD (third line of Table III). 
Thus, the previous exercise implies that the differ- 
ence between LMD and our approach is entirely due 
to the choice of mass scale, which in our case is 
determined sclf-consistcntly. Therefore, while low- 
est meson dominance might be a reasonable strat- 
egy when little information about the correlator is 
known, with the present status of the HLBL, it 
seems more justified to take full advantage of the 
available information. 



A. Extension of the DIP ansatz 

The interpolator considered in the previous Sec- 
tion provided an estimate for the HLBL that fulfilled 
the leading long and short-distance constraints from 
F 7r o 7 » 7 .. However, as we stated in the introduction, 
in this paper we want to investigate the dependence 
on the low-energy parameters (3 and 7 and estimate 
the impact of the pion away from its mass shell. 

The parameters $ and 7 were computed in Sec- 
tion [Ell Eqs. ([69)) and (|70f . In order to understand 
their potential importance for the HLBL, in Fig. 5 
we show the shape of the weight functions entering 
the dominant ajjli. As pointed out in the rel- 
evant contributions are highly peaked at low values 
of momenta, < Q 2 < 0.5 GeV. Therefore, this sug- 
gests that more information on (very) low energies 
can help improve the determination of a M . 

In order to check the stability of the results ob- 
tained so far, in particular Eqs. (|86[) . we consider a 
generalization of the DIP interpolator by adding an 
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Model 


wiGi 


W2G2 


W3G3 


W4G4 


a M 


LMD 
DIP* 
DIP mp 


+0.015 
+0.018(3) 
+0.015 


+0.042 
+0.034(4) 
+0.043 


+0.0016 
+0.0016 
+0.0016 


-0.0002 
-0.0002 
-0.0002 


7.3 • 10~ 10 
6.7(3) • 10~ 10 
7.35 • 10" 10 



TABLE III: Determination of a M with the DIP parametrization: in the second row, the mass scale is determined 
dynamically while, in the third row, my = m p . Comparison is made with the LMD model. We explicitly show the 
contribution from the different generalized form factors d. 




FIG. 5: The integrand of Eq. (78) . 



extra pole: 



K( qi ,q 2 ) = l + J2\i 

i 

2 

+E 7 * 



4 



9 9 1 9 9 

qf — mq q% — m t 

9 9 

(qf - m 2 )(<zf - mf) 



Notice that this increase in the number of poles, 
which is also natural from a large-iV c perspective, 
makes Eq. ([88)) still compatible with the general form 
of Eq. ([751) , now with 



f(q 2 ) = ! + ][>+][>, 



2 » 

q z — m~ 



m(q 2 ) 



Aj + rji 



q 2 — mf 



(89) 
(90) 



Therefore, the HLBL contribution can still be com- 
puted with the double integrals of Eqs. ([75| and ([79")) . 
By expanding Eq. (|88[) at low and high energies and 
matching it afterwards to the information coming 
from QCD short distances and low-energy (experi- 
mental and holographic) input, the constraints we 
would like to impose take the following form: 



2 X 

e4 



(a) Long-distance constraints: 

(-1.76 + 0.22) GeV~ 2 , (exp.) (91) 

(3.33 ± 0.32) GeV -4 , (pred.) (92) 

(-2.84 + 0.21) GeV" 4 , (pred.)(93) 



E^ 

* — ' mj 



2 X 



which correspond to the low-energy parame- 
ters d, and 7 defined in the previous Sec- 
tion. The value for d is taken from experiment, 
while /3 and 7 are estimated from the spread of 
values predicted by the HW1, HW2 and SW 
models and given in Eqs. (|69|) and (|70l) . 



= 



(b) Short-distance constraints 
2 2 

1 + 2 ^ + e 

i i 

2 

^m 2 (Ai + rn) 

i 

i + X> 



2 f 2 



3 



-Xo 



(94) 
(95) 
(96) 



The first two are the ones already discussed in the 
previous Sections, namely the absence of constant 
terms and the matching of the l/Q 2 coefficient. The 
last condition comes from Eq. ((731) . 

We have already mentioned that the DIP inter- 
polator (and its extensions) do not depend, by con- 
struction, on the pion momentum. While this al- 
lows us to use Eqs. ((75)) and (JTHJ), we want to note 
at this point that this simplifying assumption im- 
plies that Eq. above and the Brodsky-Lepage 
constraint cannot be satisfied simultaneously (un- 
less Xo = 0). However, there are reasons to expect 
that the Brodsky-Lepage constraint plays a minor 
role in HLBL. First and foremost, the kernel entering 
(g— 2) M is peaked at low energies (see Fig. [3]) and pre- 
sumably information on the high-energy region will 
have a negligible impact on the anomalous magnetic 
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moment. In this respect, it is quite reassuring that 
our analysis, in which the F v o~*^* low-energy region 
is approximately captured by Eqs. (f9"Tj) - (|9li|) . turns 
out to be in very good agreement with analyses that 
comply with the Brodsky-Lepage constraint [H, @ . 

The detailed solution of the system of Eqs. (pTTj) - 
(|96[) is given in Appendix B. The main conclusion of 
the analyses is that there is a remarkable stability 
when cither or 7 arc added as input. Furthermore, 
the uncertainty on (g— 2)^ is dominated by d, while 
the effect from j3 and 7 is subleading. Our estimate 
for the pion contribution to the HLBL is 

a M = 6.54(25) • 1CT 10 . (97) 



B. Impact of \o on (g — 2) M 

Since xo cannot be accessed experimentally, one 
has to resort to nonperturbative techniques to es- 
timate it. However, the different available esti- 
mates show strong discrepancies: while recent re- 
sults based on QCD sum rules [45| or exclusive B 
decays [46| seem to favor values hovering around 
Xo ~ 2 — 4 GeV -2 , estimates based on the axial 
anomaly [U EU , Pade approximants [49| or the orig- 
inal sum rule determination of (50| point at much 
higher values, xo ~ 8 — 11 GeV -2 . A recent result 
in the context of holographic QCD [fnj also pointed 
out the possibility that xo ~ 0. Most of these esti- 
mates are given without reference to the renormal- 
ization scale. However, its low-energy scale evolu- 
tion was studied in [52| and does not account for the 
discrepancies found. Therefore, the best one can do 
at present is to consider a conservative < xo < 8.9 
GcV -2 . 

The result in Eq. (|9"T)) is obtained assuming that 
Xo ~ (1 — 3) GeV -2 , which seems to be the range 
of values favored by our interpolator (see Table V 
in Appendix B). However, larger values of xo ~ 9 
GeV -2 , which at present cannot be excluded, might 
induce a shift of at most 15%. Given the precision 
presently needed for the (g — 2)^, we believe that a 
better understanding of xo is essential. 



V. CONCLUSIONS 

The muon anomalous magnetic moment is one of 
the most precisely measured quantities in particle 
physics, and thus a key parameter to test the appear- 
ance of physics beyond the Standard Model. With 
experimental accuracies expected to soon reach the 
10 -10 level, it is essential to have the theoretical un- 
certainties in the hadronic light-by-light scattering 
contribution well under control. 



In this paper we have studied the electromagnetic 
pion form factor i 7 ' 7r o 77 with a set of holographic 
models. We have shown that only a restricted set 
of the holographic models studied in Section [Til *• e 
HW1, HW2 and a modified version of SW, are able 
to predict a value of the linear slope compatible with 
the experimental data. Within this restricted set of 
models we have then obtained, in Section IIII1 pre- 
dictions for the quadratic slopes 

j3 = 3.33(32) GeV -4 , 

7 = 2.84(21) GcV -4 . (98) 

Notice that, because of the holographic models used, 
the previous predictions are sensitive to the way chi- 
ral symmetry breaking is implemented and whether 
the spectrum is Regge-like or not (SW versus HW). 
This sensitivity is contained in the quoted uncer- 
tainty. 

We have then introduced an interpolator for the 
form factor that encodes both short distances (com- 
ing from QCD) and long distances (the quadratic 
slopes /? and 7 predicted above together with the 
experimental linear slope d) and computed the pion 
exchange HLBL contribution to the (g— 2)^. Follow- 
ing 0, we have also estimated the effects due to a 
departure from the pion-polc approximation through 
a new short-distance constraint for i 7 ' 7r o* 7 » 7 », which 
is controlled by the parameter xo- Our final number 
for the (pion exchange) HLBL contribution is 

af = (6.54 ±0.25) • 10 -10 , (99) 

which depends on the slope (d) and curvature {/3 
and 7) of F 7r o 7 » 7 » as well as on the parameter xo- 
Our analysis shows that, while /3 and 7 are impor- 
tant ingredients, the uncertainty is currently domi- 
nated by the experimental accuracy on d. In view of 
this, improvements on the slope d should be a pri- 
ority. Actually, this is one of the main points in the 
scientific program for the KLOE-2 proposal [ll[ at 
DA$NE in Frascati. With high enough statistics, it 
is even feasible that the holographic predictions for 
$ and 7 could be tested. 

Another sizable source of uncertainty comes from 
the parameter xo- Our result in Eq. (|M1) assumes 
< xo is 3 GcV -2 . While numbers hovering around 
Xo ~ 2 GeV -2 seem to be favored by our analysis 
and also by different theoretical estimations, higher 
values (xo ~ 9 GeV -2 ) have also been reported in 
the literature. These higher values for xo would in- 
duce a (10 — 15)% shift on the HLBL contribution. 
Therefore, we want to emphasize the need to put 
stronger bounds on the current value of xo- In par- 
ticular, a lattice determination of xo would be ex- 
tremely useful. 
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Finally, the effects of model dependence in our 
approach are evaluated in two ways: first, through 
the spread of values for the low-energy quartic terms 
in i 7 ! 7r o 7 » 7 * , coming from the different holographic 
models; and second, using different versions for the 
parametrization of i 7 ' 7I .o 7 * 7 ». The resulting uncer- 
tainties are negligible, as it can be seen in Table IV. 

We have also compared our analysis with LMD 
models. In Section [TTT] we have investigated the issue 
of LMD in i 7 ' 7r o 7 * 7 ,» with holographic models, which 
contain a full spectrum of resonances. Our conclu- 
sion is that LMD is a sound approximation. This is 
confirmed later on in our evaluation of the HLBL: 
the mass scale of our interpolator, which is deter- 
mined dynamically, is always close to m p . In other 
words, the LMD observed in (g—2)^ is a consequence 
of LMD in the i 7 ! lr o 7 . 7 », a feature that holographic 
realizations successfully predict. A natural conse- 
quence of LMD is that the curvature 7 of the pion 
form factor can be estimated in terms of the slope a 



7 = 



2.97 GeV" 



(100) 



which agrees well with the holographic prediction. 
Also $ can be estimated from LMD, but an addi- 
tional ingredient has to be provided. If one uses 
the first short- distance constraint coming from the 
DIP interpolator and combines it with lowest vector 



dominance setting my 



one finds 



2mi 



3.10 GeV" 



(101) 



again in remarkable agreement with the holographic 
prediction. This also means that the short- distance 
constraint in the DIP interpolators is satisfied to a 
very good approximation by the low-energy values 
predicted by holographic models. This compatibility 
between long and short distances is of considerable 
importance. Since the previous results are mostly 
based on p-dominance in conjunction with our DIP 
interpolator, we conjecture that they might have a 
broader applicability, not just in ir° — > 77 but also in 
leptonic and semileptonic decays, e.g., Kl — > fJ- + ^~ 
and Kl —> Tr°e + e~ . In particular, it opens up a way 
to estimate $ in Kl — > 7r°e + e~). 



However, when it comes to the evaluation of the 
HLBL, within the current precision LMD is prob- 
ably not good enough. We have shown that the 
discrepancy between our result and the one from 
LMD models originates entirely from the different 
choices of the effective mass scale: while LMD as- 
sumes my = m p , in our approach the mass is de- 
termined from the matching conditions. Since this 
deviation is comparable to the typical uncertainty on 
HLBL, assuming my = m p does not seem justified. 
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Appendix A 

The generalized form factors Gi entering 
Eqs. (|78I79P arc given by 

Gi(x,y) = --^-/(-x 2 )F,o r7 ,(- 2/ 2 ,0) , 

N g(—x 2 ) 
G ^ y) = ~ 12^ ml ^rA-V'V* 

G 4 (x,y) = /(0)f 7r o 7 * 7 .(-x 2 ,-y 2 ) - G 3 {x,y) 



12tt 2 / w 

while the weight factors Wi take the form 



(Al) 
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wi(x,y) 



6m 2 xy(y 2 + m%) 



[4(y 2 - 2m 2 )(x 2 - y 2 ) 2 log(l + A(0, x, y)) 



4 6 ,-4 4 
i (I Jo i 6\ 24 o 4 2 , 2/ 1 S 

+ ni + y ) — x y — 3x y + 



2m% 



m- 



Amlx 2 y 2 - (x 2 ~ y 2 ) 2 m ) 



i 2/ i , i\£,i x ) 4/2 r, 2\2s(?/j . / 2 o 2\( 2 2 o 2/ 2 , 2 \ \ S \ x )£ \V) 

+ y (y -im^—^-x (y -2m*) ^— j + (y - 2m*)(x y -2m fi (x +y )) ^4 



(A2) 



U! 2 (m,x,y) 



6m 2 xy(y 2 + m 2 



[4(y 2 - 2m 2 )(x 2 - y 2 ) 2 log(l + A(0, a, y)) 



4(y 2 - 2mi )(m 4 + (x 2 - y 2 ) 2 + 2m 2 (x 2 + y 2 )) Iog(l + A(m, x, y)) 



+ ^ (m 4 + 2m 2 (, 2 + y 2 ) + (x 4 + y 4 ) - 2x 2 y 2 ) Vm - ^x 2 (y 2 - 2m 2 )£(y) + ^x 2 y 4 



(x 2 - y 2 ) 2 r] - to 6 - 3m 4 (x 2 + y 2 ) - 3m 2 (x 4 + y 4 ) - 2m 2 x 2 y 2 } 1 - 



(A3) 



w 3 (m,x,y) 



6m 2 xy 



[4(x 2 - y 2 )(m 2 (y 2 - x 2 ) + 2x 2 y 2 ) log(l + A(0, x, y)) + 4 log(l + A(m, x, y)) 



x (m 4 m 2 + (x 2 - y 2 )(m 2 (x 2 - y 2 ) - 2x 2 y 2 ) + 2m 2 (m 2 l {x 2 + y 2 ) + x 2 y 2 )) + ra 4 (x 2 + y 2 ) 



+ 2m 2 (x 4 + y 4 ) - m 2 x 2 y 2 - m 2 (m 2 + 2x 2 + y 2 )£(x) - to 2 (to 2 - 3x 2 + 2y 2 )£(y) - m 2 £(x)£(y) 
+ {(* 2 + y 2 )% - (m 2 +x 2 + y 2 )r, m } (x 2 - £(x)) + {(y 2 - 3x 2 )r /0 - (m 2 + y 2 - 3x 2 )r? m } (y 2 - £(y))] . 

(A4) 



r 



In the previous expressions we have defined 

(m 2 + x 2 + y 2 - 77 m )(x 2 - £(x))(y 2 - £(y)) 



X(m,x,y) 



8m 2 x 2 y 2 



and 



- ^ 4 + 4m 2 ; 



(A5) 



7/ m = -y/ (m 2 + x 2 + y 2 ) 2 — 4x 2 y 2 , (A6) 



(A7) 



parameters /3, 7 and \o- This we will do by vary- 
ing the combination of constraints to be imposed on 
Eq. (|Bil . 



First let us concentrate on the parameters f3. The 
system of equations we have to solve is then 



^ ml 



= -1.76(22) GeV" 2 , (exp.) 



Appendix B 



We start from the generalized interpolator: 



Qi 



t \ o 9*9 2 

v q\ — mi q% — vn~ 



9 2 



(qf - m 2 )(q 2 - mf) 



, (Bl) 



whose expansion at low and high energies leads to 
the system of constraints of Eqs. (| 9 II) - ([96]) . Our pur- 
pose is to investigate the impact of the low-energy 



* — ' rru 

i 1 
2 2 

l + 2^A l + ^r ?l 

i i 

2 

'Y^m 2 {\ i +77,) 



3.33(32) GcV~ 4 , (pred.) 



4- 2 / 2 



(B2) 



In order to ensure a solution, we will set 7712 = m p = 
0.775 GeV and A2 = 0. The interpolator, therefore, 
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reduces to 
K(i)(qi,q 2 ) 



= 1 + A 



3? 



! 1 2 
L 92 

2 2 
<7l<?2 



'{qi-mj){ql-mf) 



(B3) 



In order to check for stability, we will also consider 
the setting in which rri2 = m p = 0.775 GeV and 
772 = 0. The corresponding interpolator then reads 



K{2) (91,92) = 1 + E A '< 



95 



9? 



9 9 



9 9 

9r95 



(9? 



m 2 ) 



(B4) 



The main source of uncertainty comes from the slope 
a, while the effect of $ or 7 is subleading. Thus, an 
improvement of the CELLO data on the slope of the 
7r°77 form factor can help decrease the theoretical 
uncertainty on the HLBL. 

The previous estimate for a M most probably over- 
estimates the uncertainty Notice that the results 
for DIP^ ' give lower values of a M (due to a strong 
suppression of w\G\) and with bigger uncertainties 
from a. This might be a consequence of the strong 
correlation between a and 7 in the interpolator we 

are using and, therefore, such an uncertainty might 

(1 2) 

be misleading. Relying on DIPV ' ; and DIP„ only, 
a more realistic estimate is 

a p = 6.54(25) • 1CT 10 . 



(B7) 



The advantage of dealing with two generalized ex- 
pressions is that we will be able to assess the de- 
pendence of our results on the specific form of the 
interpolator. We will denote the previous setting as 

(i) 

DIP - , where the superscript refers to the dif- 

p;m2—m p 

ferent interpolator used (■K'(i) or K^)) and the sub- 
script emphasizes that 777,2 has been fixed and that 
$ has been used as input. The results for the HLBL 
are shown in the second and third rows of Table IV 
and can be compared with the result from DIP&, 
which corresponds to the single resonance analysis 
(see Table III) in the main text. The stability of the 
results is remarkable, as can be seen from the fact 
that the dependence on the parametrization (-f^(i) 
or K(2)) turns out to be very mild. 

An analogous analysis can be done for 7. The 
only difference is that now the system of matching 
equations looks like 



Next we turn to the evaluation of the impact of 
Xo- For this we will keep -fsT(i) and K( 2 ) but, instead 
of setting 7772 = Tn p , we will add the short-distance 



constraint involving xo to Eq. (|B2 



e4 



-1.76(22) GeV" 2 , (exp.) 



\ = 3.33(32) GeV" 4 , (pred.) 



1 + 2^^ + ^77, = 

i i 

2 

^m 2 (Aj +77,) = 

i 

i+X> = 



47T 2 / 2 

3 

47T 2 / 2 



-Xo 



(B8) 



— ' mf 

1 1 

2 A 
* — ^ in 



-2.84(21) GeV" 4 , (pred.) 



1+2 E A '+E f '« 

i i 

2 

E m ^ Ai+ ^) 



. 



47T 2 / 2 



The entries in Table IV are DIPV. m 



= 6.4(5) • 10" 



-1.76(22) GeV 2 , (exp.) We will denote the resulting interpolators as 

DIP\_ . Accordingly, the addition of the \o con- 
straint to Eqs. (|B5|) will be denoted as DIPlj . It 
is by no means trivial that those new systems of 
constraints will have a solution at all. It crucially 
depends on the interpolator. In this respect, no- 
tice that the DIP ansatz has the right kinematics to 
comply with Eq. (|73[) . which is highly nontrivial. 

The results are shown in the last three lines of 
Table IV, where we scanned the solution over a wide 
interval for xo i mostly covering all the values quoted 
in the literature so far. The range for xo is indicated 
in each case. While, in principle, we took < xa < 
8.9 GeV -2 , wherever the interval is shorter, it means 
that no solution was available beyond that point. In 
particular, DIPl 1 ^ gave no solution at all in the full 
interval. 

Taken at face value, the results of Table IV in- 
dicate a potential (10 — 15)% effect due to the our 



(B5) 



. From the 

results, we see again that the different determina- 
tions agree within errors, but somehow the uncer- 
tainty associated with a has grown. 

A conservative estimate, taking into account the 
first five rows of Table IV would give 



(B6) 
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Model 




W2G2 


W3G3 


W4G4 




Ulr s 

/3;m 2 =m p 


+0.018(3) 
+0.014 


+0.034(4) 
+0.037 


+0.0016 
+0.0015 


-0.0002 
—0.0002 


6.7(3) • 10~ 10 
6.52(15)(10) • 10 


DIP (2) 

P;m2—m p 


+0.014 


+0.037 


+0.0015 


-0.0002 


6.55(21)(6) ■ lO" 10 


DIPl 1 ' 
DIP (2) 

7;m 2 — m p 


+0.004 
+0.002 


+0.047 
+0.047 


+0.0015 
+0.0015 


-0.0002 
-0.0002 


6.09(81) (9) • 10 -10 
6.21(77)(7) • 10" 10 


DIP'. 1 ' 

,3;0<X0<8.9 

DIP (2) 

/3;0<xo<4-4 
■p)Tp(2) 


[+0.003; +0.047] 
[+0.002; +0.027] 
[+0.006; +0.035] 


[+0.043; +0.022] 
[+0.044; +0.025] 
[+0.043; +0.022] 


[+0.0015; +0.0016] 

+0.0015 
[+0.0015; +0.0016] 


-0.0002 
-0.0002 
-0.0002 


[5.9; 8.9] ■ 10" 10 
[6.0; 6.7] • lO" 10 
[6.3; 7.3] • 10" 10 



TABLE IV: Determination of a M with the generalized DIP parametrizations and their comparison with the DIP^ 
model. We explicitly show the contribution from the different generalized form factors d . The errors in parentheses 
are the ones induced by a followed by /3 or j. For the last 3 rows, where xo is inside a range, we show the values for 
the endpoints. 



Model 


a (GeV" 2 ) 


P (GeV -4 ) 


7 (GeV" 4 ) 


Xo (GeV" 2 ) 


DlPa 


-1.76* 


2.67 


4.25 


2.42 


DIP mp 


-1.35 


1.73 


2.25 


1.66 


DIPl 1 ' 

/3;m. 2 — m p 


-1.76* 


3.33* 


3.78 


1.61 


DIP (2) 

p;m2-m p 


-1.76* 


3.33* 


3.88 


1.69 


DIP* 1 ' 


-1.76* 


4.56 


2.84* 


-0.81 


DIP ( ~ 2) 

7;m 2 = m p 


-1.76* 


5.23 


2.84* 


-0.74 


DIP'. 1 ' 

/3,0<xo<8.9 


-1.76* 


3.33* 


[3.10; -5- 10 5 ] 


[0;8.9]* 


DIP (2) 

/3,0<xo<4.4 


-1.76* 


3.33* 


[3.19; -3.18] 


[0;4.4]* 


"T)Tp( 2 ) 


-1.76* 


[5.47,-18] 


2.84* 


[0; 6.4]* 



TABLE V: Predicted low-energy parameters from the different parametrizations. Asterisked quantities are input of 
the different DIP interpolators. For the last 3 rows, where xo is inside a range, we show the values for the endpoints. 
Notice from the last three lines that large values for \o are clearly disfavored. 



present knowledge of the parameter xo- The way 
to read this result is that the off-shellness of the 
pion can induce, at most, a 15% shift on values ob- 
tained within the pion-polc approximation. While 
the uncertainty on xo can only be reduced with a 
reliable calculation (probably from the lattice), we 
can still investigate which are the preferred values 
for xo for the different interpolators we have been 
dealing with. The results are collected in Table V. 
Notice that lower values for xo are preferred (first six 
rows). In contrast, large values for xo give wrong 



predictions for the low-energy parameters (see the 
last 3 rows) and hence they are clearly disfavored. 
While this might be indicative, it is certainly not 
conclusive. It is therefore essential to have a bet- 
ter understanding of xo- Incidentally, Table V also 

(1 2) 

shows that DIP^, ' ' predicts not only slightly large 
values for j3 but also negative central values for xo- 

(1 2) 

This seems to confirm that the results of DIP. 1 ; 

7 

might be unreliable. 
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